Intrinsic Friction of Monolayers Adsorbed on Solid Surfaces. 
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ABSTRACT 

We overview recent results on intrinsic frictional properties of adsorbed monolayers, 
composed of mobile hard-core particles undergoing continuous exchanges with a vapor 
phase. In terms of a dynamical master equation approach we determine the velocity 
of a biased impure molecule - the tracer particle (TP), constrained to move inside the 
adsorbed monolayer probing its frictional properties, define the frictional forces exerted by 
the monolayer on the TP, as well as the particles density distribution in the monolayer. 

INTRODUCTION 

Monolayers emerging on solid surfaces exposed to a vapor phase are important in dif- 
ferent backgrounds, including such technological and material processing operations as, 
e.g., coating, gluing or lubrication. Knowledge of their intrinsic frictional properties is 
important for understanding of different transport processes taking place within molecular 
films, film's stability, as well as spreading of ultrathin liquid films on solid surfaces [1], 
spontaneous or forced dewetting of monolayers [2,3] or island formation [4]. 

Since the early works of Langmuir, much effort has been invested in the analysis of 
the equilibrium properties of the adsorbed films [5]. Significant analytical results have 
been obtained predicting different phase transitions and ordering phenomena. As well, 
some approximate results have been obtained for both dynamics of isolated non-interacting 
adatoms on corrugated surfaces and collective diffusion, describing spreading of the macro- 
scopic density fluctuations in interacting adsorbates [6]. 

Another important aspect of dynamical behavior concerns tracer diffusion in adsorbates, 
which is observed experimentally in STM or field ion measurements and provides a useful 
information about adsorbate's viscosity or intrinsic friction. This problem is not only a 
challenging question in its own right due to emerging non-trivial, essentially cooperative 



behavior, but is also crucial for understanding of various dynamical processes taking place 
on solid surfaces. Most of available theoretical studies of tracer diffusion in adsorbed layers 
(see, e.g., [7]) exclude, however, the possibility of particles exchanges with the vapor, which 
limits their application to the analysis of behavior in realistic systems. 

Here we focus on this important problem and provide a theoretical description of the 
properties of tracer diffusion in adsorbed monolayers in contact with a vapor phase. More 
specifically, the system we consider consists of (a) a solid substrate, which is modelled in 
a usual fashion as a regular lattice of adsorption sites; (b) a monolayer of adsorbed, mo- 
bile hard-core particles in contact with a vapor and (c) a single hard-core tracer particle 
(TP). We suppose that the monolayer particles move randomly along the lattice by per- 
forming symmetric hopping motion between the neighboring lattice sites, which process 
is constrained by mutual hard-core interactions, and may desorb from and adsorb onto 
the lattice from the vapor with some prescribed rates dependent on the vapor pressure, 
temperature and the interactions with the solid substrate. In contrast, the tracer particle 
is constrained to move along the lattice only, (i.e. it can not desorb to the vapor), and 
is subject to a constant external force of an arbitrary magnitude E. Hence, the TP per- 
forms a biased random walk, constrained by the hard-core interactions with the monolayer 
particles, and always remains within the monolayer, probing its frictional properties. 

The questions we address here are the following: First, we aim to determine the force- 
velocity relation, i.e., the dependence of the TP terminal velocity V tr (E) on the magnitude 
of the applied external force. Next, we study the form of the force- velocity relation in the 
limit of a vanishingly small external bias. This allows us, in particular, to show that the 
frictional force exerted on the TP by the monolayer particles is viscous, and to evaluate 
the corresponding friction coefficient. Lastly, we analyze how the biased TP perturbs the 
particles density distribution in the monolayer; we proceed to show that there are stationary 
density profiles around the TP, which mirror a remarkable cooperative behavior. We also 
consider the case of monolayers sandwiched between two solid surfaces and demonstrate 
that here the cooperative phenomena are dramatically enhanced. Detailed account of these 
results can be found in the original papers [10-13]. 

We finally remark that our analysis can be viewed from a different perspective. On 
the one hand, the system under study represents a certain generalization of the "tracer 
diffusion in a hard-core lattice gas" problem (see, e.g., [7]) to the case where the random 
walk performed by the TP is biased and the number of particles in the monolayer is not 
explicitly conserved, due to exchanges with the reservoir. We recall that even this, by 
now classic model constitutes a many-body problem for which no exact general expression 
of the tracer diffusion coefficient D tr is known. On the other hand, our model provides 
a novel example of the so called "dynamical percolation" models (see, e.g. Ref. [8] for a 
review), in which the monolayer particles act as a fluctuating environment, which hinders 
the motion of an impure molecule, say, a charge carrier. Lastly, we note that the model 



under study can be thought of as some simplified picture of the stagnant layers emerging in 
liquids being in contact with a solid body. It is well known (see, e.g. Ref.[9]) that liquids 
in close vicinity of a solid interfaces - at distances within a few molecular diameters, 
do possess completely different physical properties compared to those characterizing the 
bulk phase. In this "stagnant" region, in which an intrinsically disordered liquid phase 
is spanned by and contends with the ordering potential of the solid, liquid's viscosity is 
drastically enhanced and transport processes are essentially hindered. Thus our model can 
be viewed as a two-level approximate model of this challenging physical system, in which 
the reservoir mimics the bulk fluid phase with very rapid transport, while the adsorbed 
monolayer represents the stagnant layer emerging on the solid-liquid interface. 

THE MODEL AND THE EVOLUTION EQUATIONS 

Consider a two-dimensional solid surface with some concentration of adsorption sites, 
which is brought in contact with a reservoir containing identic, electrically neutral particles 
- a vapor phase (Fig.5), maintained at a constant pressure. For simplicity of exposition, we 
assume here that adsorption sites form a regular square lattice of spacing a. We suppose 
next that the reservoir particles may adsorb onto any vacant adsorption site at a fixed rate 
f /t*, which rate depends on the vapor pressure and the energy gain due to the adsorption 
event. Further on, the adsorbed particles may move randomly along the lattice by hopping 
at a rate 1 /At* to any of 4 neighboring adsorption sites, which process is constrained by 
hard-core exclusion preventing multiple occupancy of any of the sites. Lastly, the adsorbed 
particles may desorb from the lattice back to the reservoir at rate g/r*, which is dependent 
on the barrier against desorption. Both / and g are site and environment independent. 

Further on, at t = we introduce at the lattice origin an extra hard-core particle, 
whose motion we would like to follow; position of this particle at time t is denoted as R tr . 
We stipulate that the TP can not desorb from the lattice and that it is subject to some 
external driving force, which favors its jumps into a preferential direction. Physically, this 
situation may be realized if the TP is charged, while the monolayer particles are neutral, 
and the whole system is subject to external electric field. The TP transition probabilities 
are determined by: 



where (3 is the reciprocal temperature, (E ■ e) stands for the scalar product, the charge of 
the TP is set equal to unity and the sum with the subscript \x denotes summation over all 
possible orientations of the vector e M ; that is, \x = {±1,±2}. We suppose here that the 
external force E is oriented according to the unit vector ei. 






Figure 1: Adsorbed monolayer in contact with a vapor. Grey spheres denote the monolayer 
(vapor) particles; the smaller black sphere stands for the driven TP. 

In the general d- dimensional case, the time evolution of P(R tr , 77; t) - the joint probabil- 
ity of finding at time t the TP at the site R tr and all adsorbed particles in the configuration 
77, where 77 = {?7(R)} denotes the entire set of the occupation variables of different sites, 
obeys the following Master equation [10-13] 

1 d 

o t p(R tr , m t) = — E { p ( Rt - ^ *) - p ( Rt - * £ ) } 

+ - - ^( R tr)) P ( R tr - e M , 77; t) ~ (1 - 77(R tr + e„)) P(R tr , V, *) } 

+ {(l-^(r))P(Rt r ,r7 r ;t)-r7(r)P(R tr ,r7;t)} 

+ {^(r)P(Rtr,7/ r ;t)-(l-7 / (r))P(R tr ,7 / ;t)}. (2) 

where 77 r,i/ is the configuration obtained from 77 by the Kawasaki-type exchange of the 
occupation variables of two neighboring sites r and r + e„, and fj v - a configuration obtained 
from the original rj by the replacement r/(r) — > 1 — 77(1-), which corresponds to the Glauber- 
type flip of the occupation variable due to the adsorption/desorption events. 

The mean velocity V tr (t) of the TP can be obtained by multiplying both sides of Eq.(2) 
by (Rtr • ei) and summing over all possible configurations (R tr ,?7), which yields 

V »® = J t £ (Rtr ' e i) P ( R tr> V, t) = ^{pi (l - fc(e i; *)) - p-! (l - fc(e_ i; t)) }, (3) 

Rtr,?7 



where 

fc(A; t)=J2 ^( R tr + A)P(Rt P , V, t) (4) 

R-tr,?? 

is the probability of having at time t an adsorbed particle at position A, defined in the 
frame of reference moving with the TP. In other words, fc(A; t) can be thought of as being 
the density profile in the adsorbed monolayer as seen from the moving TP. 

Note that Eq.(3) signifies that the TP velocity is dependent on the monolayer particles 
density in its immediate vicinity. If the monolayer is perfectly stirred, i.e., if k(X;t) = p s 
everywhere, where p s = f / (g+f) is the Langmuir adsorption isotherm [5], one would obtain 
from Eq.(3) a trivial mean-field result = (pi — — p s )(a/r), which states that 

the only effect of the medium on the TP dynamics is that its jump time r is renormalized 
by (1 — Ps) 1 - We proceed to show, however, that k(X;t) is different from p s everywhere, 
except for |A| — > oo. This means that the TP strongly perturbs the monolayer. 

The method of solution of Eq.(2) and (3) has been amply discussed in [10-13]. Here 
we merely present the results. 



RESULTS FOR ONE-DIMENSIONAL MONOLAYERS 



We focus here on the one-dimensional lattice, which is appropriate to adsorption on 
polymer chains [14], and on the limit t — > oo. In this case, the stationary particles density 
profile, as seen from steadily moving TP, has the following form: 



k n = k(X) = p s + K±exp(^ — cr\n\/X± j , X — an, n e Z, 
where the characteristic lengths A± obey 



A ± = =F a In 1 



A 1 + A_, + 2(f + g) T J (A, + A_, + 2(f + <?)) - 4A t A_ 
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the amplitudes K± are given by 

A 1 - A_! 



K 4 



Ps 



l A_ 1 -A 1 exp(-a/X + ) , Ps A^ exp(-<r/A_) - A 1 ' 

the terminal velocity V tr = <j{Ai — A_i)/2r*, while A 1 and A_i obey: 



A x -A. 



and 



A 1 = 1 + 
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Pit 
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Figure 2: The TP velocity as a function of the adsorption probability / at different values of 
g. The TP hopping probabilities are p\ = 0.6 and p-\ = 0.4. The solid lines give the analytical 
solution while the squares denote the Monte-Carlo simulations results. Upper curves correspond 
to g = 0.8, the intermediate - to g = 0.5 and the lower - to g = 0.3, respectively. 



Note that A_ > A + , and consequently, the local density past the TP approaches its non- 
perturbed value p s slower than in front of it. Next, K + is always positive, while K_ < 0; 
this means that the density profile is a non-monotonous function of A and is characterized 
by a jammed region in front of the TP, in which the local density is higher than p s , and a 
depleted region past the TP in which the density is lower than p s . 

For arbitrary values of p, f and g the parameters A ±1 , Eqs.(8) and (9), and consequently, 
V tr can be determined only numerically (see Figs. 2 to 4). However, V tr can be found 
analytically in the limit E — > 0. In the leading in E order, we find 



tr 



(10) 



which relation is an analog of the Stokes formula for driven motion in a ID adsorbed 
monolayer undergoing continuous particles exchanges with the vapor phase and signifies 
that the frictional force exerted on the TP by the monolayer particles is viscous. The 
friction coefficient ( is given explicitly by 



c 



2r 



1 + 



PsT* 



/?a 2 (l - p s ) L r(f + g)l+ ^1 + 2(1 + r*(l - Ps )/r)/(f + g) 



(11) 



Note that the friction coefficient in Eq.(ll) can be written down as the sum of two contri- 
butions ( = ( cm + (coop- The first one, ( cm = 2t/ f3a 2 (l — p s ) is a typical mean-field result 



and corresponds to a perfectly homogeneous monolayer. The second one, 



C = ?rV. 1 (12) 

Wop (5o\\ - p a ){f + g) 1 + yiT2(T+ r*(l - Ps )/r)/(/ + ^ ' 1 J 

has a more complicated origin and reflects a cooperative behavior emerging in the mono- 
layer, associated with the formation of inhomogeneous density profiles (see Fig. 3) - the 
formation of a "traffic jam" in front of the TP and a "depleted" region past the TP (for 
more details, see [10]). The characteristic lengths of these two regions as well as the am- 
plitudes K± depend on the magnitude of the TP velocity; on the other hand, the TP 
velocity is itself dependent on the density profiles, in virtue of Eq.(3). This results in an 
intricate interplay between the jamming effect of the TP and smoothening of the created 
inhomogeneities by diffusive processes. Note also that cooperative behavior becomes most 
prominent in the conserved particle number limit. Setting /, g — > 0, while keeping their 
ratio fixed (which insures that p s stays constant), one notices that Ccoop gets infinitely large. 
As a matter of fact, in such a situation no stationary density profiles around the TP exist; 
the size of both the "traffic jam" and depleted regions grow in proportion to the TP mean 
displacement X tr (t) ~ \ft [15]. 
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Figure 3: Density profile around stationary moving TP for / = 0.1, g = 0.3 and p = 0.98. The 
solid line is the plot of the analytical solution. Squares are the MC results. 

Consider finally the situation with E — 0, in which case the terminal velocity vanishes 
and one expects conventional diffusive motion such that 



Xl(t) = 2 Art, 



(13) 



where D tr is the diffusion coefficient. We can evaluate D tr if we assume that the Einstein 



tr 



relation D tr — f3/( holds, which yields 
cr\l- Ps 



D 



tr 



1 + 



PsT 



- 'I' ' r(f + g)l + ^1 + 2(1 + r *(l-p s )/r)/(/ + g) 
Monte Carlo simulations evidently confirm our prediction for D tr (see Fig. 4). 
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Figure 4: The TP self-diffusion coefficient vs the adsorption probability /. Notations and values 
of g are the same as in Figs. 2. 



RESULTS FOR TWO-DIMENSIONAL ADSORBED MONOLAYERS 

We turn now to the case of a 2D substrate. In this case, the general solution for the 
density profiles, in the frame of reference moving with the TP, reads [10-13]: 

k(X) = fc m , na = Ps + ^2 A " h M V -" F ^ - p s {A 1 - A_i)(Vi - V_i)F nii „ 2 | (15) 

V 

with 

(A \ nx l 2 f°° / \ 

F nuna = \ -=±\ J e-*I ni (2a' 1 v / A 1 A_ 1 tj I n2 (2a- 1 A 2 t) dt, (16) 

where l n (z) stands for the modified Bessel function, the coefficients A v are determined 
implicitly as the solution of the following system of three non-linear matrix equations 



Vi/ = {±1, 2}, A v = 1 + — Pv \l-ps- p.(Ai ~ 

t \ det C 



(17) 



/A 1 V- 1 F ei - a A_iViF ei A 2 V_ 2 F ei 

c=l .i,v ./•'« , .1 V;/;, , - n A 2 v_ 2 F e _ 1 ). (is) 

V AiV^iFe, A^ViF^ A 2 V_ 2 F e2 - a, 

the matrix C v stands for the matrix obtained from C by replacing the v-t\i column by the 
column-vector F, 

(Vi - V_!)F ei \ 

(Vi-V-OFe., , (19) 

(Vi - V_!)F e2 / 

while h(e u ) are expressed in terms of A v as 

h(e v ) = (l-p 3 ) + -^-(l-A v ). (20) 

Lastly, we find that the TP velocity obeys 

a r 4r* pi det (7i - p_i det (7_i -i - 1 



a r 4r 

Vtr = -(pi 1 + Ps — 

r It 



(21) 



detC > 

Non-linear Eqs.(15) are quite complex and their explicit solution can not be obtained 
analytically. Typical density profiles are depicted in Fig. 5 and are characterized by a 
condensed, traffic-jam-like region in front of the stationary moving TP, and a depleted by 
particles region past the TP. 

The asymptotical behavior of the density profiles at large distances from the TP is quite 
spectacular. In front of the TP, the deviation h n>0 = k n $ — p s always decays exponentially 
with the distance: 

expf-n/A+J x r f— — 2 ] 

„■/» • A + = '" (xiji-^ + Vd-- 42 ) - A ' A -j\ 

(22) 

On contrary, the behavior past the TP depends qualitatively on the physical situation. 
In the general case when exchanges with the vapor phase are allowed, the decay of the 
density profiles is still exponential: 



h_ n>0 ~ F_ 



exp 



(23) 

In the conserved particles number limit, which can be realized for the monolayers 
sandwiched in a narrow gap between two solid surfaces, A_ becomes infinitely large and 
the deviation of the particle density from the equilibrium value p s follows 

h -^-M 1 + ln + °{^- (24) 



Figure 5: Typical density profiles in 2D case. The arrow shows the direction of external field. 



Remarkably enough, in this case the correlations between the TP position and the particles 
distribution vanish algebraically slow with the distance! This implies that mixing of the 
monolayer is not efficient enough to prevent the appearance of the quasi-long-range order 
and the medium "remembers" the passage of the TP on a long time and space scale. 
Now, in the limit f3E C 1 we find a Stokes-type formula of the form V tr ~ E/(, where 

c _ 4r jl | T * ^ 1 (25) 

/Ml -Ps)y t (j + g + 1 + T *( 1 _ Ps y T ^ - x) J ' 

with 

(26) 

Note that we again are able to single out two physically meaningful contributions to the 
friction coefficient (. Namely, the first term on the rhs of Eq.(25) is just the mean- field- 
type result corresponding to a perfectly stirred monolayer, in which correlations between 
the TP and the monolayer particles are discarded. The second term, similarly to the ID 
case, mirrors the cooperative behavior emerging in the monolayer and is associated with 
the backflow effects. In contrast to the ID case, however, the contribution stemming out 
of the cooperative effects remains finite in the conserved particles limit. 



Lastly, we estimate the TP diffusion coefficient D tr = (3 1 ( 1 as 



D tr = -7-{l -p 8 )<l 



4r v ,h ' \ r A(f + g+l + T*(l-p s )/r)C(x)-l + (3p s -l)T*/T) 

(27) 

Note that setting / and g equal to zero, while assuming that f /g = p s / (1 — p s ), we recover 
from the last equation the classical result due to Nakazato and Kitahara (see [7]), which 
is exact in the limits p s <ti 1 and p s ~ 1, and serves as a very good approximation for the 
self-diffusion coefficient in hard-core lattice gases of arbitrary density [7]. 



CONCLUSIONS 



To conclude, we have studied analytically the intrinsic frictional properties of adsorbed 
monolayers, composed of mobile hard-core particles undergoing continuous exchanges with 
the vapor. We have derived a system of coupled equations describing evolution of the 
density profiles in the adsorbed monolayer, as seen from the moving tracer, and its velocity 
V tr . We have shown that the density profile around the TP is strongly inhomogeneous: the 
local density of the adsorbed particles in front of the TP is higher than the average and 
approaches the average value as an exponential function of the distance from the TP. On the 
other hand, past the TP the local density is always lower than the average, and depending 
on whether the number of particles is explicitly conserved or not, the local density past 
the TP may tend to the average value either as an exponential or even as an algebraic 
function of the distance. The latter reveals especially strong memory effects and strong 
correlations between the particle distribution in the environment and the TP position. 
Next, we have derived a general force- velocity relation, which defines the TP terminal 
velocity for arbitrary applied fields and arbitrary values of other system parameters. We 
have demonstrated next that in the limit of a vanishingly small external bias this relation 
attains a simple, but physically meaningful form of the Stokes formula, which signifies that 
in this limit the frictional force exerted on the TP by the monolayer particles is viscous. 
Corresponding friction coefficient has been also explicitly determined. In addition, we 
estimated the self-diffusion coefficient of the tracer in the absence of the field. 
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